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Model  for  the  Computation 
us  Army  Corps  of  Time-Steadv  Nearshore  Currents 

of  Engineers^ 

by  Bradley  Johnson 


PURPOSE:  The  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  herein  introduces 
a  new  model  for  the  computation  of  depth-dependent  time-steady  currents  in  the  surf  zone. 

BACKGROUND:  Under  the  assumptions  of  longshore  uniformity,  monochromatic  and 
unidirectional  waves,  one-dimensional  (1-D)  models  have  predicted  the  steady  depth- integrated 
longshore  current  with  reasonable  accuracy  in  the  regions  of  intense  breaking  (e.g.,  Bowen  1969; 
Longuet-Higgins  1970;  Kraus  and  Earsen  1991).  The  addition  of  random  wave  fields  with 
directional  spread  (Battjes  1974;  Thornton  and  Guza  1986)  produce  similar  results  but  indicate  a 
difference  in  the  magnitude  of  the  horizontal  mixing  terms.  It  is,  however,  necessary  to  also  predict 
the  vertical  structure  of  time-steady  currents  in  the  surf  zone  to  reliably  compute  sediment  transport 
rates.  Indeed,  concentrations  of  sediment  are  found  to  be  orders  of  magnitude  larger  in  the  near-bed 
boundary  layer  where  fluid  velocities  are,  in  general,  smaller  than  the  depth-averaged  value.  It  is, 
therefore,  of  practical  importance  to  predict  the  vertical  variation  of  nearshore  currents. 

The  phase  resolving  equations  of  fluid  motion  may  be  best  suited  for  a  physically  based  prediction  of 
sediment  transport.  Ein  and  Liu  (1998),  for  instance,  solved  the  Reynolds  equations  with  an 
algebraic  nonlinear  Reynolds  stress  model.  Karambas  and  Koutitas  (2002)  and  Kobayashi  and 
Johnson  (2001)  used  Boussinesq  and  nonlinear  shallow-water  formulations,  respectively.  These 
time-dependent  numerical  models,  however,  are  not  suitable  to  most  practical  engineering  problems 
due  to  the  prohibitively  large  computation  time. 

The  new  model  NEARHYDS  efficiently  predicts  wave  height,  setup,  bed  ripple  formation, 
boundary  layer  thickness,  and  depth-dependent  steady  cross-shore  and  longshore  currents.  The 
solution  of  the  time-averaged  equations  of  energy  balance  and  alongshore  momentum  is 
implemented  with  modular  coding  that  is  flexible  and  well  suited  for  the  iterative  process  required  to 
include  the  effects  of  interaction.  NEARHYDS  is  a  compiled  collection  of  Matlab  routines  and  a 
graphical  user  interface  complete  with  an  on-line  user’s  guide  and  is  suitable  for  the  computation  of 
steady  nearshore  hydrodynamics  on  beaches  with  essentially  straight  and  parallel  contours.  The 
compiled  version  runs  on  a  Windows  PC,  and  the  hydrodynamic  predictions  for  a  typical  surf  zone 
domain  are  completed  in  several  seconds.  This  physics-based  two-dimensional  (2-D)  model  has  been 
funded  through  the  Advanced  Nearshore  Circulation  work  unit  of  the  U.S.  Army  Corps  of  Engineers 
Navigation  Systems  Research  Program  and  will  form  the  basis  for  a  more  general  version  that  allows 
alongshore  variability. 

No  mechanism  exists  within  the  formulation  of  the  equations  for  the  generation  of  low-frequency 
standing  or  progressive  waves  in  the  cross-shore.  These  low  frequency  motions,  however,  may  play 
an  important  role  in  cross-shore  profile  development.  While  the  wave-current  interaction  is  included 
in  the  bottom  boundary  layer  and  stress  analysis,  the  effect  of  the  currents  on  the  short  waves  was 
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not  considered.  As  justifieation,  it  was  found  that  the  deviation  from  Snell’s  law  due  to  longshore 
eurrent  is  exeeedingly  small,  and  the  eross-shore  eurrents  are  eonsidered  to  be  small  relative  to  the 
wave  phase  speed  on  alongshore  uniform  eoasts. 

MODEL  DESCRIPTION:  The  deeay  in  wave  height  of  obliquely  ineident  waves  is  the  primary 
generation  meehanism  for  longshore  eurrents;  undertow  is  generated  through  the  mass  flux  due  to 
the  presenee  of  waves,  a  surfaee  roller,  and  the  depth-inhomogeneity  of  the  radiation  stress.  Wave 
breaking  is  ineorporated  into  the  momentum  equations  through  the  solution  of  a  suitable  phase- 
averaged  energy  equation  in  the  eross-shore  direetion.  With  an  estimate  of  the  energy  flux,  the 
radiation  stresses  due  to  the  organized  wave  motions  ean  be  quantified  and  ineluded  as  foreing  in  the 
time-averaged  and  depth-averaged  momentum  equations.  The  gradients  of  the  radiation  stress  in  the 
eross-shore  direetion  are  eounteraeted,  primarily,  by  a  slope  in  the  mean  free  surfaee  position, 
termed  the  setup.  The  alongshore  momentum  equation  is  a  balanee  of  the  radiation  stress  gradient 
with  momentum  mixing  and  bottom  shear  stress.  The  depth-variation  of  the  steady  eurrents  are  based 
on  horizontal  time-averaged  Reynolds  equations  where  use  is  made  of  the  solutions  of  the  energy 
and  depth-integrated  momentum  equations  as  well  as  physieally-based  elosures.  First,  the  solution  of 
the  phase-averaged  energy  equation  for  the  waves  is  presented,  where  two  approaehes  are  provided: 
a  root-mean-square  (RMS)  representation  and  a  wave-by-wave  (WBW)  representation.  Comparisons 
of  wave  height  are  presented  as  verifieation  of  these  methods.  Seeond,  the  solution  of  the  eross-shore 
momentum  balanee  is  deseribed,  and  results  in  the  predietions  of  setup  and  steady  eross-shore 
eurrents.  Finally,  the  longshore  momentum  balanee  is  presented  that  provides  the  longshore  eurrent. 
Comparisons  with  lab  data  are  shown  as  a  gauge  of  the  model  aeeuraey,  and  more  thorough  vali¬ 
dation  with  field  data  will  be  eompleted  as  the  development  proeess  proeeeds. 

Phase-Averaged  Energy  Equation.  The  energy  equation,  time-averaged  over  many  shortwave 
periods,  may  be  expressed  as 


^E,=-D,  (1) 

ax 

in  whieh  x  is  the  eross-shore  eoordinate  positive  onshore;  Ep  =  time-averaged  energy  flux  per  unit 

width;  and  Dp  =  time-averaged  energy  dissipation  rate  due  to  wave  breaking,  whieh  needs  to  be 

estimated  empirieally  in  this  model.  The  model  NEARHYDS  solves  Equation  1  with  either  the 
RMS  wave  height  representation  presented  by  Battjes  and  Jansen  (1978)  and  ealibrated  by  Battjes 
and  Stive  (1985)  or  expressed  as  the  superposition  of  regular  waves  (WBW)  (e.g..  Dally  1992). 

RMS  representation.  The  time-averaged  energy  flux  per  unit  width  for  random  waves  ean  be 
represented  with  the  RMS  wave  height  and  peak  period  as 

—  1  7 

EF  =  ^PsH^,nC^cosa  (2) 

where  pg  is  the  unit  weight  of  water;  =  ^/8  cr  with  a  =  standard  deviation  of  the  free-surfaee 

elevation;  n  is  the  ratio  of  the  wave  group  speed  to  linear  phase  speed;  =  phase  veloeity  based 
on  Tp ;  and  a  =  speetral  peak  angle  (Figure  1). 
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Figure  1.  Variable  definition  sketch 


The  wave  breaking  dissipation  formula  by  Battjes  and  Stive  (1985)  is  given  by 

D',  =  \psf,QHl  (3) 

where  is  the  spectral  peak  frequency  ;  Q  is  the  local  fraction  of  breaking  waves  governed  by  a 
Rayleigh  wave  height  distribution  truncated  with  the  maximum  wave  height  determined  using 
the  Miche  criteria. 


The  energy  equation  (1)  is  integrated  numerically  using  Equations  2  and  3,  resulting  in  an  estimate 
of  the  cross-shore  distribution  of  wave  height.  Recalling  that  the  energy  equation  is  the  first  step  in 
the  solution  of  the  time-averaged  momentum  balance,  the  radiation  stresses  are  then  expressed  using 
linear  theory  as 


S„=\pgH^ 


2 

rms 


n(cos^  a+l)-  — 


n  sin  a  cos  a 


(4) 
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Linear  wave  superposition  method.  Alternatively,  the  energy  flux  ean  be  expressed  as  the 
superposition  of  monochromatic  waves  under  the  assumption  of  narrow  handedness  in  frequency. 
Individual  wave  contributions  are  governed  by 


-  1  2 

EFi  =  -pgHi  nCp  cosa 

O 


i  =  \:N 


(5) 


where  H.  =  wave  height  for  the  wave,  and  N  is  the  total  number  of  waves.  The  values  of  each 
wave  height  realization  at  the  offshore  boundary  are  given  in  terms  of  the  exceedance  probability 


r 

r  0 

H.=H  1 

(in 

i  rms  .4 

(6) 


where  P^.  is  a  single  realization  of  the  random  variable  P^  that  is  uniformly  distributed  between  zero 

and  one.  The  distribution  of  individual  wave  heights  at  the  seaward  boundary  approaches  a  Rayleigh 
distribution  as  A  ^  oo 


In  contrast  to  the  RMS  representation,  the  breaking  dissipation  is  zero  for  a  single  wave  realization 
until  the  Miche  incipient  breaking  criterion  is  met.  Thereafter,  the  dissipation  is  assumed  to  be  in 
proportion  to  the  difference  between  the  wave  energy  flux  and  an  empirical  stable  energy  flux,  Ep  , 
associated  with  a  stable  wave  height  of  OAd 

(7) 


where  d  is  the  still- water  depth. 


The  energy  equation  (1)  is  integrated  numerically  for  each  offshore  wave  condition  with  energy 
fluxes  Equation  5  and  breaking  dissipation  Equation  7  resulting  in  an  estimate  of  the  cross-shore 
distribution  for  each  individual  wave  height.  The  excess  momentum  fluxes  are  then  expressed  as 
aggregates  of  the  individual  cross-shore  wave  height  distributions; 


c  1  ^ 


n(cos^  a+l)-  — 


TV  2  1  ^ 

— 'y  pgnsmacosa — V  7/f 


(8) 


Predicted  wave  height  compared  to  iaboratory  data.  Recent  laboratory  data  of  random 
waves  are  valuable  for  model  verification.  The  recirculating  Eongshore  Sediment  Transport  Eacility 
(ESTE)(Hamilton  and  Ebersole  2001)  was  designed  for  the  collection  of  detailed  free  surface  and 
velocity  data.  A  total  of  1 10  horizontal  measurement  positions  were  set  up  in  the  ESTE  basin  with 
dimensions  of  approximately  30  m  by  25  m  in  plan.  The  analysis  of  the  hydrodynamic  and 
bathymetric  conditions  over  this  region  demonstrates  that  the  modeling  of  the  tank  can  justifiably  be 
simplified  by  assuming  alongshore  uniformity.  The  following  analysis  focuses  on  a  single 
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cross-shore  transect  located  near  the  center  of  the  tank.  Spilling  waves  with  - 1 8cm  and  peak 
period  T^=\ .5s  are  used  as  a  verification  case  herein. 

The  energy  equation  (1)  is  solved  across  the  domain  using  either  the  RMS  or  WBW  formulation. 
Random  wave  fields  are  characterized  by  smoothly  varying  wave  heights  due  to  the  range  of 
breakpoints,  and  a  first  order  explicit  method  is  used  to  numerically  integrate  Equation  1 .  This 
simple  procedure  is  considered  to  be  accurate  enough  in  light  of  the  empirical  parameterizations  and 
the  smooth  variation  of  energy  across  the  surf  zone.  Figure  2  shows  the  cross-shore  variation  of 

for  the  two  formulations  as  well  as  laboratory  data.  The  data  collection  was  repeated  a  total  of  10 
times  resulting  in  10  measurements  of  wave  height  at  each  of  the  10  cross-shore  measurement 
locations.  All  of  the  data  are  included  in  Figure  2  to  demonstrate  the  degree  of  scatter.  For  each  of 
the  numerical  solutions,  empirical  parameters  have  been  chosen  based  on  accepted  recom¬ 
mendations,  and  no  attempt  has  been  made  to  calibrate  the  models  to  fit  data.  The  difference 
between  the  solutions  by  these  two  formulations  is  due  to  the  unequal  breaking  dissipation,  where 
the  RMS  method  predicts  a  larger  dissipation  in  the  deeper  region  of  the  domain.  The  RMS  method 
underpredicted  the  wave  height  in  the  outer  surf  zone,  but  agreed  well  with  data  in  the  mid-surf  On 
the  other  hand,  the  WBW  method  matches  well  with  data  in  the  offshore  region  of  the  data  but 
overpredicts  wave  height  in  the  mid-  and  inner-surf  zones. 


Figure  2.  Bathymetry  and  comparisons  of  model  wave  height  predictions  to  measured  data 
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Cross-Shore  Momentum  Equation.  The  time-averaged  cross-shore  momentum  equation  for  a 
long  straight  coast  is  given  by 


dx 


—2 


P  dz- 


Pg 
2  dx 


-pgh 


dx  dx 


dz-T^^ 


(9) 


where  u  is  the  cross-shore  velocity;  P  is  total  pressure;  h  is  the  mean  water  depth;  rj  is  the  mean 

free-surface  position;  is  a  Reynolds  stress;  and  r®  is  the  bottom  shear  stress.  The  left  hand  side 

of  Equation  9  is  made  up  of  the  forcing  terms  that  arise  from  the  organized  wave  motion  as  well  as 
current  contributions.  The  time-steady  currents  may  play  a  role  in  momentum  mixing  in  the  surf 
zone;  it  is  the  case,  however,  that  the  rigorous  inclusion  of  these  effects  is  complicated  in  the  trough 
to  crest  regions  and  requires  assumptions  about  the  wave  shapes  and  velocity  profiles.  It  is  therefore 
assumed  that  the  currents  do  not  play  a  vital  role,  and  the  velocities  in  Equation  9  are  assumed  to  be 
wave  induced.  Eurthermore,  it  has  been  shown  that  the  mean  bed  shear  stress  and  the  turbulent 
mixing  are  small  relative  to  the  forcing  and  free-surface  slope  (e.g.,  Eonguet-Higgins  and  Stewart 
1964).  Incorporating  these  simplifications,  then,  along  with  a  linear  representation  of  the  wave 
induced  velocities  and  pressures  yields  the  familiar  cross-shore  balance; 

^  =  (10) 

dx  dx 

The  momentum  equation  (10)  is  then  solved  numerically  for  the  cross-shore  variation  of  the  setup, 
77,  given  the  known  radiation  stress  from  Equations  4  or  8.  A  first-order  explicit  scheme  is  con¬ 
sidered  to  be  accurate  enough  for  this  slowly  varying  function  of  cross-shore  distance. 

To  predict  the  depth  variation  of  the  steady  cross-shore  velocity,  use  is  made  of  an  eddy  viscosity 
closure.  The  depth-dependent  time-averaged  cross-shore  Reynolds  equation  is  integrated  twice 
vertically  to  arrive  at  the  variation,  and  boundary  conditions  are  employed  to  determine  the  inte¬ 
gration  constants  and  close  the  problem.  It  is  convenient  to  decompose  the  total  velocities  into  steady 
and  time-varying  components: 


u  =  u  +  u  ;  v  =  v-l-v  ;  w  =  w+w 


(11) 


such  that  all  time-varying  portions  have  a  diminishing  time  average  i.e.,  w  =  v  =  iv  =  0  .  It  is  assumed, 
in  the  following,  that  the  steady  vertical  current  w  is  negligible  relative  to  the  horizontal  currents. 
Also,  the  cross-shore  gradients  of  integrated  turbulent  shear  are  assumed  to  be  small  (Stive  and  Wind 
1982),  and  the  specification  of  the  vertical  gradient  of  the  wave  stress  uw  is  given  according  to 
Rivero  and  Arcilla  (1995).  The  equation  that  governs  the  current  profile,  then,  is  given  as 


dz 


(w 


du 


1  d  —2  —2.  d  -2  drj_ 

-^(w  )  +  —  u 

2  dx  dx  dx 


(12) 


6 


ERDC/CHL  CHETN-VI-39 
December  2003 


The  simplest  solution  of  the  general  equation  (12)  is  the  speeifieation  of  a  depth- invariant  turbulent 
vertical  eddy  viscosity  .  In  fact,  Garcez  Faria  et  al.  (2000)  concluded  that  a  depth-dependent  eddy 
viscosity  did  not  improve  predictions  of  the  undertow  when  compared  with  field  data.  In  the  model 

NEARHYDS,  the  eddy  viscosity  is  expressed  as  .  If  the  forcing  /„  is  also  con¬ 

sidered  to  be  constant  over  depth,  the  required  two  integrations  yield  a  quadratic  profde; 


u  = 


/ 

c,z 

+  —  +  C2 


(13) 


where  z'  -  z  +  d  is  a  new  local  coordinate  system  with  z'  -0  at  the  bottom,  and  Cj  and  Cj  can  be 

specified  through  two  boundary  conditions.  Alternatively,  the  constants  can  be  determined  with  one 
boundary  condition  and  a  statement  of  mass  conservation,  which  is  the  closure  employed  in 

NEARHYDS. 


The  first  constant  Cj  is  expressed  through  a  balance  of  the  bottom  shear  stress  and  the  eddy  viscosity 

closure  with  u  .  resulting  in  Cj  =  Tjp .  The  average  bottom  shear  stress  is  solved  with  a 

boundary  layer  analysis  that  is  beyond  the  scope  of  this  publication,  but  combines  aspects  of  Grant 
and  Madsen  (1979)  with  ripple  prediction  of  Nielsen  (1992)  and  wave  friction  factors  outlined  by 
Soulsby  (1997).  The  end  result  is  the  determination  of  a  boundary  layer  thickness  S  and  a  shear 
stress  applied  to  the  interior  portion  of  the  fluid.  The  boundary  layer  solution  must  be  solved  itera¬ 
tively,  but  it  converges  within  about  three  iterations. 

It  is  common  to  determine  the  second  constant  with  a  statement  of  mass  conservation  whereby  a 
finite  mass  flux  due  to  the  presence  of  waves  is  balanced  by  a  steady  seaward  directed  current  under 
the  trough  level.  Kennedy  et  al.  (1998)  approximated  the  volume  flux  through  the  time  averaged 
continuity  equation  and  linear  long  wave  resulting  in  an  expression  for  the  average  undertow  below 
trough  level; 


Ur 


coscr 


(14) 


The  final  unknown  can  then  be  determined  by  equating  the  mass  flux  and  the  assumed  form  of  u 
expressed  in  Equation  13  integrated  from  the  bed  to  the  trough  level,  resulting  in 


r  —  d 

‘-2  -  Ur  ^ 


6k 


(15) 


where  d,  is  the  distance  from  the  bed  to  the  trough  level,  approximated  as  A  -  .  Substitution  of 

the  constants  into  Equation  13  yields 


u  —  1J+- 


]_d^ 

2v/^  dx 


W  —2  ~ 

~  +  Ur  +  gV 


(16) 
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where  u  is  approximated  as  in  the  formulation  of  . 


The  solution  of  the  undertow  problem  requires  the  speeifieation  of  the  so-ealled  radiation  stress 

terms  if'  and  in  Equation  16.  The  simplest  solution  is  the  use  of  depth  invariant  linear  long-wave 
foreing 


g  ^  ^  g  , 

u  = ,  ;  w  =  ,  = — 

\  h  \  k  dx 


One  of  the  stress  terms,  for  example,  is  then  expressed  as 


-2_S^2_gH, 


2 

rms 


(17) 


(18) 


where  -  VScr  has  been  used. 

Thus  the  depth  variation  of  the  steady  eross-shore  flow  is  solved  with  the  use  of  a  loeal  mass 
conservation,  a  shear  stress  boundary  condition,  and  linear  long-wave  theory  to  quantify  radiation 
stresses. 

Alongshore  Momentum  Equation.  The  time-averaged  momentum  equation  for  a  long  straight 
coast  is  a  balance  of  the  organized  wave  forcing  and  current  momentum  with  turbulent  mixing  and 
bottom  shear  stress. 


puvdz  =  ^j\  dz-Ty  (19) 

dx-’^t 

where  v  is  the  longshore  current.  The  left-hand-side  of  Equation  19  represents  the  combined  effects 
of  wave  forcing  and  current  dispersion.  As  in  the  cross-shore  momentum  equation,  the  steady 
currents  may  play  a  role  in  the  distribution  of  momentum  in  the  surf  zone.  These  effects  are 
neglected,  however,  due  to  the  added  complexity  and  ambiguities  associated  with  a  more  detailed 
analysis.  The  assumption,  therefore,  is  that  the  current  effects  are  small,  and  that  the  wave-induced 
forcing  dominates  this  term  and  can  be  quantified  with  linear  theory: 


puv  dz  □  puv  dz  -  S 

j  2a  JZu  ^ 


(20) 


where  is  the  radiation  stress  given  by  Equations  4  or  8. 

A  common  Reynolds-type  expression  is  employed  in  the  specification  of  the  mixing  term  and  it  is 
assumed  that  the  mixing  stress  can  be  expressed  in  terms  of  the  near-bottom  velocity  gradient.  The 
eddy  viscosity  is  typically  assumed  to  be  depth  invariant  for  the  solution  of  the  alongshore 
momentum  problem,  and  then  the  mixing  shear  term  in  Equation  19  is  expressed: 


8 


ERDC/CHL  CHETN-VI-39 
December  2003 


dx 


rri  d 

f  ,  dVg] 

II 

1.^ 

phv, — - 
dx 

\  ox  y 

(21) 


The  eddy  viscosity  is  expressed  as  in  Longuet-Higgins  (1970),  |  |  where  N  is  an 

empirical  dimensionless  coefficient,  and  |  \  is  the  distance  to  the  shoreline. 


While  the  bottom  shear  stress  could  justifiably  be  neglected  in  the  depth  integrated  cross-shore 
momentum  balance,  the  near  bed  stress  is  the  primary  force  counteracting  the  radiation  stress 
gradient  in  the  longshore  balance.  Use  is  made  of  the  quadratic  friction  law  to  define  the  time- 
averaged  stress 


<  =  pc j  \  Ug  I 


(22) 


where  the  subscript  S  indicates  that  the  velocities  are  taken  near  the  bottom  but  outside  of  the 
boundary  layer.  The  time-averaging  in  Equations  22  is  complicated  by  the  necessary  detailed 
knowledge  of  magnitudes  and  phasing  for  both  the  steady  and  time-varying  velocities.  Conse¬ 
quently,  the  square  wave  approximation  of  Nishimura  (1982)  is  employed  to  give  a  simplified 
expression  for  the  shear  stress.  In  the  interest  of  brevity,  the  equations  are  not  presented. 


Utilizing  the  previous  simplifications.  Equation  19  reduces  to 


dx 


dx 


f 


phv, 


V 


dx  , 


(23) 


The  simplified  alongshore  momentum  equation  thus  reduces  to  a  second-order  differential  equation. 
Central  differences  are  used  to  solve  for  the  near-bed  longshore  velocity  Vg  with  known  radiation 

stresses.  A  central  difference  representation  of  the  mixing  term  necessitates  the  solution  of  Equa¬ 
tion  23  as  a  system  of  implicit  linear  finite-difference  equations.  The  offshore  boundary  condition  is 
specified  by  the  user,  and  the  shoreline  boundary  condition  is  developed  by  solving  Equation  23 
without  momentum  mixing.  Note  that  the  undertow  solution  must  be  completed  prior  to  an  along¬ 
shore  momentum  balance  due  to  the  inclusion  of  the  cross-shore  velocity  in  the  bottom  shear  stress 
formulation. 


As  in  the  cross-shore  flow,  prediction  of  the  depth  variation  of  the  steady  longshore  velocity  is  made 
with  the  use  of  an  eddy  viscosity  closure.  The  depth-dependent  time-averaged  alongshore  Reynolds 
equation  is  integrated  twice  vertically  to  determine  the  steady  current  profile.  Boundary  conditions 
are  then  employed  to  determine  the  constants  of  integration.  Stive  and  Wind  (1982)  demonstrated 
that  the  horizontal  turbulent  stress  term  was  small,  and  thus  it  is  neglected.  Also,  the  vertical  wave- 
induced  flows  are  assumed  to  be  phase  shifted  relative  to  the  horizontal  flows  such  that  vw  is  small. 
Neglecting  the  steady  term  that  includes  the  steady  vertical  current  is  justifiable  for  small  bottom 
slopes,  and  the  simplified  equation  that  governs  the  current  profile  is  expressed 
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dz  dz 


+  VU 


(24) 


Where  V  is  the  depth-averaged  longshore  veloeity  under  the  trough  level.  Again,  the  right-hand- 
side  of  Equation  24  is  made  up  of  eurrent  and  wave  forcing  components  and  is  assumed  to  be 
independent  of  depth.  The  implicit  nature  of  this  equation  dictates  the  iterative  solution  employed  in 
NEARHYDS,  Proceeding  as  in  the  cross-shore  problem,  Equation  24  is  integrated  twice,  generating 
two  constants  of  integration.  As  previously,  the  eddy  viscosity  is  assumed  to  be  depth-uniform  and 
identical  to  the  formulation  used  in  the  undertow  problem.  The  first  constant  is  specified  with  the 
match  of  turbulent  stress  and  the  bottom  shear  stress  obtained  from  the  solution  of  the  depth- 
integrated  alongshore  momentum  equation  (23).  The  second  constant  is  determined  from  the 
computed  near-bottom  velocity  resulting  in 


-  -  1  a 


v  =  v^  + 


2v,,  dx 


[ur^  +  VU^Z'^  + 


PK 


(25) 


where  the  short-wave  forcing  term  vu  is  quantified  according  to  linear  long-wave  theory. 

PREDICTED  NEARSHORE  CURRENTS  COMPARED  TO  LAB  DATA:  The  ESTF  test  case 
of  spilling  waves  that  was  used  for  the  comparison  of  predicted  wave  heights  to  data  is  also  used  to 
gauge  the  accuracy  of  the  computed  setup  as  well  as  predicted  current  magnitudes  and  depth 
profiles.  The  waves  with  an  offshore  wave  angle  of  10  deg  generated  a  longshore  current  and  a 
seaward  directed  return  current  over  the  movable  bed.  A  time  series  of  velocities  was  collected  at 
10  positions  throughout  the  water  column  at  each  cross-shore  position  as  shown  in  Figure  3.  The 
portion  of  the  column  covered  by  the  dense  vertical  spacing  of  measurement  locations  is  used  to 
verify  the  cross-shore  and  longshore  current  models.  The  nearshore  currents  were  computed  using 
standard  values  for  the  coefficient  of  friction  =0.01  and  the  horizontal  mixing  coefficient 

N  -  0.04 .  No  attempt  has  been  made  to  achieve  a  better  fit  with  data  by  tailoring  these  factors. 


Figure  3.  Vertical  measuring  locations  Test  1H,  case  3 
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The  time-average  surface  position  is  shown  in  Figure  4.  As  in  the  presentation  of  the  wave  height 
data,  all  of  the  measurements  are  shown  at  each  position  in  order  to  demonstrate  the  degree  of 
scatter.  The  differences  in  the  RMS  and  WBW  methods  are  consistent  with  the  previously  shown 
wave  height  computations  considering  that  mean  free-surface  slope  is  a  result  of  the  decay  in  . 

It  is  expected,  then,  that  the  RMS  formulation  will  exhibit  a  greater  setup  through  most  of  the 
domain;  the  greatest  slope  in  the  mean  free-surface  for  the  WBW  method,  on  the  other  hand,  will  be 
in  the  nearshore  where  the  wave  heights  drop  off  quickly. 


Figure  4.  Mean  free-surface  position  data  and  predictions 

The  model  predictions  for  the  magnitude  and  profile  of  the  undertow  are  compared  with  data 
measured  at  1 0  cross-shore  locations  in  Figure  5 .  The  problems  of  accurate  undertow  prediction  are 
compounded  by  a  reliance  on  the  total  mass  flux  parameterization,  which,  in  turn,  is  dependent  on 
the  solution  of  the  energy  equation.  Despite  this  challenge,  the  undertow  magnitude  predictions  of 
the  model  compare  well  with  the  measured  data.  The  shape  of  the  current  profile  is  also  predicted 
well  with  differences  between  the  WBW  and  RMS  methods  due  to  the  contrasting  cross-shore 
distributions  of  wave  height  and  associated  radiation  stresses. 

Computed  and  measured  steady  longshore  currents  are  shown  in  Figure  6.  The  model  uses  the 
measured  (and  negative)  longshore  velocity  at  the  first  measurement  location  as  a  seaward  boundary 
condition.  It  is  interesting  to  note,  however,  that  the  specification  of  a  non-zero  boundary  condition 
on  the  seaward  boundary  has  only  a  localized  effect.  The  only  mechanism  for  non-local  generation 
and  dissipation  of  longshore  current  is  due  to  the  momentum  mixing,  and  the  influence  of  the 
seaward  boundary  drops  off  with  a  weak  mixing  effect.  The  profile  is  well  predicted  with 
NEARHYDS  and  the  difference  between  the  WBW  and  RMS  formulations  is  small. 

SUM  MARY :  The  model  NEARHYDS  solves  the  phase-averaged  energy  equation  to  determine  the 
cross-shore  distribution  of  RMS  wave  heights  and  the  associated  radiation  stresses.  In  turn,  the  wave 
stresses  act  as  a  forcing  mechanism  in  the  depth-integrated  cross-shore  and  alongshore  momentum 
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Model/Data  Undertow  Comparision 


Figure  5.  Undertow  data  and  predictions 

equations  that  include  mixing  and  bottom  friction.  Finally,  with  use  of  the  previously  described 
analyses,  the  depth  variation  of  the  currents  can  be  determined  based  on  the  horizontal  Reynolds 
equations. 

The  prediction  of  wave  height  is  shown  to  differ  based  on  the  choice  of  RMS  or  WB  W  formulations, 
although  both  methods  predict  reasonably  well.  The  differences  in  the  computed  setup  are 

consistent  with  the  varied  wave  height  predictions.  The  longshore  current  computations  agree  well 
with  data  both  in  magnitude  and  profde  shape.  The  undertow  predictions,  on  the  other  hand,  show 
greater  error  when  compared  with  data.  Additionally,  the  undertow  profile  is  shown  to  be  sensitive  to 
the  solution  of  the  energy  equation. 

Model  input  requirements  include  1-D  (cross-shore)  bathymetry  and  incident  wave  parameters 
7),,  O').  Options  include  RMS  or  WBW  formulations,  bottom  friction  formulation,  mixing 

coefficients,  as  well  as  other  solution  parameters.  NEARHYDS  computes  the  cross-shore  distri¬ 
bution  of  wave  height  and  setup.  Additionally,  the  steady  cross-shore  and  depth  dependent  current 
fields  are  computed.  Model  output  is  displayed  graphically  and  is  recorded  to  ASCII  files.  The 
graphical  user  interface  is  intended  to  make  NEARHYDS  intuitive  and  simple  to  use. 
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Model/Data  Longshore  Current  Comparision 


Figure  6.  Longshore  current  data  and  predictions 

POINT  OF  CONTACT:  For  additional  information,  or  to  obtain  a  copy  of  the  model,  eontaet 
Dr.  Bradley  Johnson,  Coastal  Proeesses  Braneh,  Coastal  and  Hydraulies  Laboratory,  U.S.  Army 
Engineer  Research  and  Development  Center,  Vicksburg,  MS  39180.  Voice:  601-634-4612,  FAX: 
601-634-4314  email:  Bradley.D.Johnson@erde.usaee.army.mil. 
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